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ABSTRACT 

Population-Ill (Pop-Ill) starformation (SF) is thought to be quenched when the mctallicity of the star- 
forming gas has reached a critical level. At high redshift, when the general intergalactic medium (IGM) 
was enriched with metals, the fraction of primordial gas already collapsed in minihalos (above the Jeans 
mass but below the mass corresponding to the efficient atomic cooling threshold for SF) was significantly 
larger than the fraction of primordial gas that had already been involved in Pop-III SF. We argue that 
this reservoir of minihalo gas remained largely in a metal-free state, until these minihalos merged into 
large systems (above the hydrogen cooling threshold) and formed stars. As a result, the era of Pop-III 
■ SF was significantly prolonged, leading to a total integrated Pop-III SF that was an order of magnitude 

^ ' larger than expected for an abrupt transition redshift. We find that the contribution of Pop-III SF to 

the reionization of hydrogen could have been significant until z ~ 10 and may have extended to redshifts 
£C) • as low as z ~ 6. Our modeling allows for gradual enrichment of the IGM, feedback from photo-ionization 

and screening of reionization by minihalos. Nevertheless, the extended epoch of Pop-III SF may result in 
complex reionization histories, where multi peaks are possible over some regions of parameter space. The 
relative contribution of Pop-III stars to reionization can be quantified and will be tested by three-year 
WMAP results: (1) if Pop-III stars do not contribute to reionization, r cs < 0.05 — 0.06 and a rapid 
reionization at z ~ 6 are expected with the mean neutral fraction quickly exceeding 50% at z ~ 8; (2) if 
the product of star formation efficiency and escape fraction for Pop-III stars is significantly larger than 
that for Pop-II stars, then a maximum r os = 0.21 is achievable; (3) in a perhaps more plausible scenario 
where the product of star formation efficiency and escape fraction for Pop-III stars is comparable to 
\Q \ that for Pop-II stars, r os = 0.09 — 0.12 would be observed, with reionization histories characterized by 

an extended ionization plateau from z — 7 — 12 where the mean neutral fraction stays in a narrow range 
of 0.1 — 0.3. This result holds regardless of the redshift where the IGM becomes enriched with metals. 
Subject headings: Cosmology: theory - Pop-III stars - Early universe - Intergalactic Medium 
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Following cosmological recombination at a redshift z ~ 10 3 , the baryonic gas filling up the universe became predom- 
C$ • inantly neutral. Given that this gas is known to be mostly ionized today, one arrives at one of the major questions in 
| current extragalactic astronomy, namely, when was the cosmic hydrogen re-ionized, and what were the sources responsi- 
• *h . ble? Recent observations have begun to provide preliminary answers to this question. The absorption spectra of SDSS 
p\ ' quasars at z ~ 6 indicate that the neutral fraction of hydrogen increases significantly at z > 6 (Becker et al. 2001; Fan et 
H i al. 2002; Cen & McDonald 2002; Pentericci et al. 2002; White et al. 2003; Wyithe & Loeb 2004; Mesinger & Haiman 2004; 
. - - « Fan et al. 2005). On the other hand, analysis of the first year data from the WMAP satellite suggests a large optical 
depth to electron scattering r cs = 0.17 ± 0.04, implying that the universe was significantly reionized at redshifts as high 
as ~ 17 ± 5 (Kogut et al. 2003; Spergel et al. 2003). While this result is not surprising in view of earlier calculations of 
the reionization history by metal- free stars (e.g. Wyithe & Loeb 2003a; Cen 2003a), it offers an intriguing empirical path 
for answering questions regarding the sources of reionization. 

Several authors have pointed out that the combination of these results can only be reconciled with an extended or multi- 
peaked reionization history (Wyithe & Loeb 2003a, c; Cen 2003a, b; Chiu, Fan, & Ostriker 2003; Haiman & Holder 2003; 
Somerville et al. 2003; Gnedin 2004). A common feature of these models is that the first generation of stars plays a 
significant role. The first stars in the universe (Pop-III) formed out of metal-free gas, relic from the big bang. Simulations 
of metal-free star formation indicate that the first stars may have high masses (M > 100M Q ), since gas cooling by 
molecular hydrogen (H2) cannot lower the gas temperature below ~ 200K (Bromm, Coppi, & Larson 2002; Abel, Bryan 
& Norman 2002). Massive, metal- free stars shine at their Eddington luminosity, Le oc M, and have roughly constant 
effective (surface) temperatures of ~ 10 5 K and lifetimes of ~ 3 x 10 6 yr, independent of their mass (Bromm, Kudritzki, 
& Loeb 2001). This implies that the number of ionizing photons produced per baryon incorporated into these stars 
(~ 8 x 10 4 ) was independent of their mass function, and larger by more than an order of magnitude than that of 
the observed Population-II (Pop-II) metal-rich stars. As a result, the emissivity of ionizing photons may be double 
peaked, leaving open the possibility of an early reionization by Pop-III stars, with an epoch of recombination following 
enrichment of the IGM above the critical threshold, before a final reionization by Pop-II stars near z ~ 6. Indeed, the 
observed enrichment of the IGM is consistent with a substantial contribution from Pop-III stars (Oh, Nollett, Madau & 
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Wasserburg 2001). On the other hand, it has been suggested that the non-monotonic behavior of these models is due 
to the assumption of a sharp transition in metallicity (Furlanetto & Loeb 2005). These authors find that the various 
feedback mechanisms lead to an extended but monotonic reionization history. 

In this paper we develop a model which follows the densities of primordial and enriched gas in collapsed systems, and 
compute co-dependent SF, reionization and metal enrichment histories. While our model is more detailed than many 
previous semi-analytic approaches in several respects, the primary qualitative difference is that Pop-Ill SF is allowed to 
proceed in gas that had already collapsed in minihalos prior to enrichment of the IGM, but which has not formed stars, 
even when the general IGM has been enriched to a high metallicity level. We motivate the physics of this conclusion, in 
addition to modeling its consequences for the reionization history. We show that the additional post enrichment Pop-Ill 
SF boosts the total integrated Pop-Ill SF by an order of magnitude in typical cases. 

Throughout the paper we assume cosmological parameters obtained through fits to WMAP data (Spergel et al. 2003). 
These include density parameter values of Q m = 0.27 in matter, = 0.044 in baryons, £!a = 0.73 in a cosmological 
constant, and a Hubble constant of Hq = 71 kms _1 Mpc _1 (or equivalently h = 0.71). For calculations of the Press- 
Schcchtcr (1974) mass function (with the modification of Shcth & Tormcn 1999) we assume a primordial power-spectrum 
with a power-law index n = 1 and the fitting formula to the exact transfer function of Cold Dark Matter, given by Bardeen 
et al. (1986). We adopt an rms amplitude of cr 8 = 0.84 for mass density fluctuations in a sphere of radius 8/i _1 Mpc. 



We make the distinction between the ionizing radiation field due to a possible early population of zero-metallicity stars, 
and the metal enriched stars observed at lower redshifts. It is thought that the primordial initial mass function favored 
massive stars (Bromm, Copi & Larson 1999, 2001; Abel, Bryan & Norman 2000; Mackey, Bromm, Hernquist 2002). The 
possible existence of this population is very important for reionization because the spectrum of these stars would result 
in an order of magnitude more ionizing photons per baryon incorporated into stars (Bromm, Kudritzki, & Loeb 2001). 
The formation of the very massive stars is suppressed as the material out of which stars form is enriched with metals. 
The fraction of the ionizing photons produced by metal-free stars depends on several unknown parameters, including the 
mixing efficiency of metals, the environments in which new stars form, and most importantly, the threshold metallicity 
above which star formation is transformed from being dominated by massive stars to a Scalo (1998) initial mass function 
(IMF). The threshold metallicity is believed to be small; Bromm ct al. (2001) argue for a threshold Zt ^'° Bh < 10~ 3 of the 
solar metallicity value. The efficiency of mixing of metal enriched outflows from star forming galaxies to the surrounding 
IGM is even more uncertain; Scannapieco, Ferrara & Madau (2002) and Ricotti & Ostriker (2004) find that the mass 
weighted mean metallicity can reach values greater than 10~ 3 of the solar value at redshifts as high as 20, and note that 
the average metallicity scales with star formation efficiency, supernova rate, and the fraction of supernova energy that is 
channeled into outflows. In a representative region of the IGM, the average metallicity increases roughly exponentially 
with redshift as it is modulated by the exponential growth in the collapse fraction of baryons at high redshifts. We 
therefore define a transition redshift z tran below which metal- rich stars with a Scalo (1998) IMF dominate the production 
rate of ionizing photons. However, because enrichment will not occur at the same time in all parts of the IGM, we allow 
for a spread in the transition redshift. We assume a Gaussian probability function with variance A2; tran , so that the 
fraction of IGM (/t ran ) that is enriched at redshift z is given by 



The effectiveness of a stellar population in ionizing hydrogen in the IGM can be parameterized in terms of the number 
of ionizing photons produced per baryon incorporated into stars (N). Throughout the paper we assume Pop-Ill stars to 
be massive (> 100M©) zero-metallicity stars with the generic spectrum calculated by Bromm, Kudritzki & Loeb (2001). 
The resulting value is A^ pr i mor d = 84168. We assume metal enriched stars (l/20th solar metallicity) to have a Scalo (1998) 
mass-function, and use spectral information from the stellar population model of Leitherer et al. (1999) 3 . This results in 
a value of 7V en rich ~ 4270. 



Throughout the paper we assume that star-formation is initiated through atomic hydrogen cooling, and so implicitly 
assume the rapid destruction of molecular hydrogen in the IGM by a back-ground of UV photons; molecular hydrogen 
had been destroyed long before the universe, on average, became significantly ionized (an ionized fraction of ~ 10 -3 ) 
(Haiman, Rees & Loeb 1997; Oh & Haiman 2003). Therefore, at the time when the IGM started to become enriched to 
the critical level for Pop-II SF there was significant gas already collapsed in minihalos (halos above the cosmological Jeans 
mass, but below the minimum mass for SF) that had virial temperatures which were too low to initiate SF (see § 4.4). If 
this gas was preserved in its primordial state until SF was initiated through formation of a sufficiently massive system, 
then the epoch of Pop-Ill SF could have been increased to redshifts significantly below z tran . 



Model spectra of star-forming galaxies were obtained from http://www.stsci.edu/science/starburst99/. 



2. IONIZING RADIATION FROM POP-III AND POP-II STARS 




(1) 
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3.1. Enrichment of Virialized Primordial Gas In Minihalos 

It is often assumed that metal enrichment of the intergalactic medium (IGM) implies metal enrichment of all gas 
in the universe (e.g., Wyithe & Loeb 2003a, b; Cen 2003a, b Furlanetto & Loeb 2005). We shall point out that this 
simple assumption is incorrect, which, as we will show, may lead to qualitatively different conclusions with regard to 
contributions of Pop-Ill stars to cosmological rcionization. A critical point to note is that metal enrichment process 
facilitated by Shockwaves carrying metals generated in stars is a strong function of the density of the gas that is being 
enriched. In particular, we will make the distinction between metal enrichment of the average IGM and metal enrichment 
of gas that is inside already formed halos. 

Let us consider metal enrichment of gas in minihalos when a blastwave of velocity U sweeps through the IGM. This 
is equivalent to the case of a gas cloud moving through a background medium with velocity U. The gas cloud may 
be subject to a host of hydrodynamical instabilities, including Kelvin-Helmholtz and Richtmeyer-Meshkov instabilities, 
through which metal mixing occurs. Several authors (Murray et al. 1993, M93 hereafter; Klein, McKee, & Colella 1994; 
Dinge 1997; Miniati et al. 1997) have shown that a non-self-gravitating gas cloud moving at the sound speed of the 
background medium gets disrupted on a time scale of the dynamic time of the cloud. Self-gravity, in our case by the 
dominant dark matter in a minihalo, could stabilize the gas cloud. M93 studied the case of self-gravitating, non-cooling 
gas clouds moving through a background medium in the context of a two-phase medium. The gas clouds in our minihalos 
are physically identical to the case considered by M93 in that there is little cooling (due to lack of cooling agents at the 
low temperature), and the gas cloud is confined by gravitational potential well of the dark matter minihalo. Therefore, we 
may draw directly upon the calculations by M93 to quantify the whether or not gas clouds in minihalos survive enrichment 
of the IGM. M93 show both analytically and by simulations that the stability condition for a gas cloud moving through 
a background medium requires the surface gravity of a cloud to be greater than 

2ttC/ 2 

where D is the density ratio of the gas cloud to the background gas, R c \ is the radius of the gas cloud. If we define 77 as 
then we find that 

„( r ) = 8 5( U r 2 ( — ) 2 /3 ( l±£ )( ^L)-4.7 (4) 

V( ' l 10km/s J { 10 (i h- 1 M Q ) { 15 A Mj { ' 

where we have assumed that the density slope near the virial radius is —2.4 (Navarro, Frenk, & White 1997); Mh is the 
mass of the minihalo within its virial radius; M r is the mass with radius r, z is redshift. M93 show at rj — 1 only 2% and 
11% of the gas is lost after 3.2 and 10 times the dynamic time of the cloud, which is about the Hubble time. However, 
they also note that even at rj = 0.25, the gas mass loss is still relatively small. 

While U might be large in the immediate vicinity of a shockwave-producing galaxy, one does not expect U to be large 
at large distances. Mori, Ferrara, & Madau (2002) show, in simulations of propagation of supernova blastwaves from 
10 s Ii~ 1 Mq galaxies at z = 9, that after more than a hundred million years the relative filling factor for regions being 
swept by shocks of velocities larger than U = (10, 30, 100)km/s is roughly (100%, 35%, 10%). We expect the velocities to 
be still smaller at the higher rcdshifts of concern here, due to enhanced cooling and larger Hubble velocity. We see from 
equation (4) that if U < 29km/s 100% of the gas within virialized regions of halos more massive than 1O 6 /i _1 M will 
have T) > 1 and be relatively un-affected by the Shockwaves. This fraction is reduced to 60% and 30% for U = lOOkm/s 
and U = 500km/s, respectively. For lower mass mini-halos (Mh ~ 10 4 5 Mq) close to the Jeans mass (§ 4.1), we find r\ ~ 1 
for a smaller velocity of U ~ lOkm/s. However n(r) is a high power of the fraction of gas in the halo that is subject to 
instabilities. Thus even for halos of mass Mh ~ 10 4 ' 5 M Q , and winds of velocity U ~ 30km/s and U ~ lOOkm/s we find 
values of -j^ ~ 0.6 and ^ ~ 0.35 yield rj ~ 1. Combining these results, we expect that most of the gas already virialized 
with minihalos will be largely unaffected by metal-carrying blastwaves and remain metal-free until its first episode of 
starformation. 

In addition to direct enrichment, there is the possibility that virialized primordial gas in mini-halos could be enriched 
during the process of mergers with other mini-halos containing enriched gas. We do not expect this to be the case for 
the reasons outlined below. Following the transition redshift, if the blastwaves that carry metals sweep through the IGM 
with a velocity of order of ~ 10 km/s, the IGM may be collisionally ionized. The collisionally heated IGM subsequently 
cools via the Compton cooling process (by the cosmic microwave background), and at the same time recombines in the 
absence ionizing sources. We find that collisionally ionized gas would be Compton-cooled and saturate at a temperature 
of 400 — 500K (depending on the initial temperature) after about a Hubble time at z — 15; at 1/4 of the Hubble time 
T > 600 K. Hence, IGM that is heated at the redshift of interest will not be able to cool below <~ 400 K by radiative 
processes. Subsequent cooling is primarily due to adiabatic expansion, which could achieve a temperature reduction by a 
factor of a few for the redshift range of interest. For a IGM of T ~ 400 K at z = 15 the Jeans mass is ~ 4 x 10 6 M Q , which 
is close to the threshold mass for atomic cooling (see equation 7 below). Therefore, metal-enriched and heated IGM is 
likely to be prevented from further accreting onto all minihalos, new and old. Thus, gas in collapsed minihalos will not 
be contaminated by enriched IGM gas and will remain metal-free, and newly formed minihalos will be gas-free. 

To be extremely conservative, let us consider the case where the metal-enriched and heated IGM is assumed to have 
cooled to the temperature of the cosmic microwave background. Gas from the IGM can enter mini-halos with M > 
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10 5 — 10 6 Mq either through mergers of smaller collapsed objects, or through accretion from the IGM. In the former case, 
the cold gas content is primordial because accretion from the IGM at the temperature of the cosmic microwave background 
is suppressed in low mass mini-halos. In the latter case we argue that newly accreted gas into more massive mini-halos is 
unlikely to be as dense as the gas already at the virial radius. This may be seen as follows. An estimate of the non-linear 
overdensity of baryons inside a virialized object is 

, Pb 1 A 6T vir \2 9 T vir 

where pb is the density of baryons inside the virialized object of temperature T v ; r , and pb and T are the background 
baryon density and temperature respectively (Barkana & Loeb 2001). At rcdshifts below z ~ 130, and in the absence of a 

heating source the temperature of an adiabatically cooling IGM is T g ~ 6.8 (^p) K (Peebles 1993). This temperature is 
substantially lower than the temperature of the cosmic microwave background, Xcmb ~ 54.6 (-^p). Thus primordial gas, 
which has not been subjected to heating accretes into a dark-matter halo with an initial temperature that is substantially 
lower than does primordial gas which has been heated by the enriching blast-wave before being Compton cooled to the 
temperature of the cosmic microwave background. From equation (5) we therefore find that the ratio of overdensities 

of primordial to enriched gas following virialization is <~ 8 (^p) • This result suggests that in the absence of cooling, 
enriched gas would be prevented from significantly penetrating the minihalo, so that mixing is expected to be inefficient. 
As a result primordial gas is not expected to become contaminated during hierarchical growth of mini-halos. 

Upon merging to form a system large enough to initiate SF, does the primordial gas in mini-halos form stars before 
becoming enriched by other SF in the galaxy? This is expected to happen to some extent. But "other" SF that occurred 
earlier would have been Pop-Ill SF. Therefore, the question is not whether Pop-Ill SF will occur or not, but rather what 
the efficiency of Pop-Ill formation will be. We do not have a quantitative answer to this at present. It seems plausible 
to expect that the first SF would occur in the largest mass clump in the galaxy. If the largest mass clump makes up a 
significant fraction of the total mass of the galaxy, then the overall SF formation efficiency may be largely determined by 
the SF efficiency within that largest mass clump. We will therefore assume, following conventional wisdom, that Pop-Ill 
SF proceeds with efficiency r\ = 0.1. 

In summary, the above discussion leads us to explore the effect of the following hypothesis on the Pop-Ill SF history, 
and on the reionization history of cosmic hydrogen, i) We assume that virialized primordial gas within mini-halos is not 
enriched by super-galactic winds, ii) Virialized primordial gas is not enriched following accretion of enriched gas from the 
IGM, or during the mergers of mini-halos. Hi) As a result, primordial virialized gas remains in its primordial state until 
such time as its host halo grows to be massive enough to initiate a burst of Pop-Ill SF. 



4. EVOLUTION OF THE DENSITY OF GAS IN HALOS 

In order to compute the effect of IGM enrichment on the SF history, we need to follow the density of enriched and 
primordial gas through the hierarchical merging of halos. This section introduces a method for following this evolution. 
During hierarchical galaxy formation, an average property of the galaxy population, such as the density of galaxies evolves 
in redshift due to newly collapsing halos, mergers of halos and accretion. We begin by discussing the critical mass-scales 
in galaxy formation. We then discuss calculation of the rate of newly collapsing halos, before describing how this quantity, 
along with the merger rates of galaxies can be used as source and sink terms in differential equations that describe the 
average evolution of the gas content of the galaxy population. 



4.1. Critical Mass Scales 

There are four critical mass scales affecting galaxy formation. The first is the Jeans mass, which corresponds to the 
smallest mass halo into which gas can accrete from the IGM. By considering the response of a baryonic overdensity to 
the potential well of a dark-matter halo one can find the mass scale at which the baryonic overdensity reaches 100. We 
take this mass scale to represent the Jeans mass (Barkana & Loeb 2001) 

The second mass scale corresponds to the virial temperature above which gas that has accreted into a halo can cool 
efficiently. We assume this scale to be determined by atomic cooling, and so to correspond to a virial temperature below 
Tmin ~2x 10 4 K. The corresponding mass is therefore 

m {l.98xW*Kj Up/ U, V 10 / ' [ ' 

where Sl z m = [1 + (f2 A /Q m )(l + z)" 3 ]" 1 , A c = 18tt 2 + 82d - 39d 2 and d = VL z m - 1 (see Barkana & Loeb 2001 for more 
details). The third mass scale becomes important following the reionization of a region of IGM at which time it is heated 
to ~ 10 4 K. The Jeans mass is increased by several orders of magnitude, and numerical simulations find that gas infall is 
suppressed in halos with T v ; r < 2.5 x 10 5 . There is some disagreement as to the exact value of the halo circular velocity 
below which gas infall is completely suppressed (e.g. Thoul & Weinberg 1996; Kitayama & Ikeuchi 2000; Quinn, Katz & 
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Efstathiou 1996, Weinberg, Hcrnquist & Katz 1997; Navarro & Steinmetz 1997; Dijkstra, Haiman, Rees & Weinberg 2004). 
In what follows we assume a post-reionization minimum temperature of T re i on ~ 2.5 x 10 5 , and note that the assumption 
of a smaller temperature will not qualitatively effect our conclusions. Thus we have 
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T • \ : 
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2.5 x 10 5 K J 
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(8) 
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Finally, the fourth characteristic mass is the non-linear mass-scale (M n i), which corresponds to the co- moving length 
scale i? n i = ( 4^°' ) 1 ^ 3 over which typical fluctuations (i.e. 1-sigma fluctuations) in the linear density field have a value 

of a(M n i) — S CT it/D(z). Here p m is the co-moving density of matter, <5 cr it is the linearly extrapolated overdensity at the 
time of collapse for a dark matter halo, and D(z) is the growth factor at redshift z. For the cosmology employed in this 
paper, and at rcdshifts z > 1, a good fit to the dependence of the non- linear mass scale on redshift is given by 



M nl (z) = 10 13 - 58 (1 + z 



^1.39 



exp 



0.18 



0.79 



(9) 



From equation (9) we see that M n i moves through the Jeans mass at z ~ 8, and the mass corresponding to the hydrogen 
cooling threshold (M m i„) at z ~ 4. The non-linear mass-scale corresponds to the typical mass of newly forming galaxies. 
Equation (9) therefore encapsulates the reason why reionization is expected to occur at redshifts between z <~ 10 and 
z ~ 5. Reionization is a self- limiting process because the non-linear mass-scale does not reach M rc i on until z ~ 2. 

4.2. Primordial Gas In Collapsed Systems At z tran 

Prior to the enrichment of the IGM with metals produced by the first stars the gas content of collapsed systems below 
the critical threshold for SF (minihalos) was metal-free. Systems above the critical threshold had an initial burst of 
metal-free SF, but then presumably became enriched internally so that subsequent SF was metal enriched. Following 
enrichment of the IGM to a mctallicity above a critical value, gas that was newly accreted from the IGM was enriched, 
and so resulted in Pop-II SF. However, the metal-free gas in minihalos would not have been enriched and would have 
moved through mergers into larger systems. When this metal-free gas became part of a galaxy above the SF threshold, 
further Pop-Ill SF would have been possible, even significantly after the IGM became enriched. Indeed, if one calculates 
the collapsed fraction of gas in galaxies at z = 15 with masses above the Jeans mass, to the collapsed fraction above the 
mass corresponding to the minimum for SF (M m ; n ) one finds that the collapsed fraction differs by a factor of 20! Figure 1 
shows the ratio of collapsed fractions above M m j n and the Jeans mass as a function of redshift. The majority of metal-free 
gas that has accreted into collapsed systems at the time of enrichment has not yet formed stars. As a result the time 
of enrichment may not correspond to the time when Pop-Ill SF ends. Indeed the fraction of the Pop-Ill SF during the 
history of the universe that has taken place by the transition redshift could be as low as the ratio shown in Figure 1. We 
therefore expect most Pop-Ill SF to take place at redshifts below z tra n- 




O20 15 10 

z 

Fig. 1. — The ratio of collapsed fractions of gas in galaxies with masses above the Jeans mass, and above M m i n respectively. 



4.3. The Rate of Newly Collapsing Dark-Matter Halos 
To follow the evolution of primordial gas in a hierarchical cosmology we first need to compute the collapse rate of new 



halos. Let -j^f(z) be the Press-Schechter (1974) mass function (number of halos with mass between M and M + dM 
per co- moving Mpc 3 ) of dark-matter halos at redshift z. The rate of change of the density of dark matter halos between 

M and M + dM is therefore dMi P d B z ■ Let us also define (Lacey & Cole 1993) dA] ^ d s t \ M the number of mergers per unit 
time of halos of mass AM with halos of mass M (forming new halos of mass Mi = M + AM) at redshift z. We find 
the component of the change in density with redshift of halos having masses between M and M + dM that is not due 
to mergers. This is given by the overall change in the density of halos with mass M, minus the change in density due 
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to merger of halos to form a new halo of mass M, plus the change in density of halos of mass M that merge with other 
halos. The collapse rate of new halos may be written 



d 2 n co \ _ d 2 n ps f 2 , A , , d 2 N, 
dMxdz ~ dMidz 



[ 2 dAM d N ™z * dn^ r 

J dAMdt Mi _ AM dz d{M 1 - AM) J dAMdt 



dAM 



d N, 



dt dn 



ps 



Mi dz dM\ 



(10) 



,2 

The value of d ^° dz equals the net change in density due to newly collapsed halos of mass M. 



4.4. The Evolution Of The Primordial Gas Content Of Galaxies 

Our next step is to find the evolution of the density of metal- free gas per halo mass dp pi i mm( i/dM. Note that this 
density refers to a density of gas that has not been in a galaxy at a time when it underwent a burst of SF. We assume 
that if a halo forms with a mass Mi larger than the Jeans mass Mj, or if the mass of a halo grows through a merger to 
become larger than Mj, then enough gas accretes inside the virial radius so that the mass of gas inside the halo becomes 
fbM, where /b is the fraction of the mass-density in baryons. Therefore, if the halo mass is larger than Mj, then the 
total density of gas inside halos in a small mass-range between Mi and Mi + 5M is 

dPi~d SM = f M *g.sM. (11) 
dMi J dMi v ; 

Our aim is to compute the effect of metal enrichment on the density of cold, metal free gas inside galaxies of mass Mi . 

Prior to enrichment of the IGM, the gas in galaxies where SF has not occurred remains metal-free. Following enrichment, 
newly accreted gas is metal enriched. In neutral regions the gas inside halos below M min is not photo-evaporated, and 
so, at least initially some halos with masses below M m i n will contain metal-free gas. This gas will move into larger halos 
during subsequent mergers. We can calculate the density of metal-free gas that remains in these low mass halos following 
enrichment of the IGM. If two halos M and AM merge to form a new larger halo Mi with a mass smaller than M m j n , 
then the average mass of metal-free gas in the new halo will equal the sum of the average masses of metal-free gas in each 
of the initial halos + %^/^). 

We can therefore write down the evolution of the density of primordial gas which we denote dPp J^°' d (recall that this 
density refers to a density of gas that has not been in a galaxy at a time when it underwent a burst of SF) 

r M i 

Mi 

+ ^ dKiidi M ,_ A „ S«-4M]' 

The first equation of (12) is valid in neutral regions, while once the region become ionized any remaining primordial gas is 
photo-evaporated out of the minihalos on the halo dynamical timescalc (~ O.IH^ 1 ), where it then becomes enriched, so 

d 2 

that 4MTdz d = 0" ^ rs * ^ ne con t ams two terms that relate to regions of IGM that have not yet been enriched. The 
first term corresponds to a source of collapsing halos into which metal-free gas accretes. The second term corresponds 
to a source of gas accumulation in halos above the Jeans mass. This accumulation is the sum of gas that was already in 
halos if the progenitor was above the Jeans mass, and gas newly accreted from the IGM if it was not. The third term 
(second line) corresponds to regions with enriched IGM and describes the movement of primordial gas from small halos 
to large halos during mergers. The fourth term describes the loss of gas density in halos of mass M that results from 
the merger of those halos to form larger systems. The primordial gas density in halos smaller than the Jeans mass or 
larger than M m ; n is zero. The upper left panels of figures 2 and 3 show the evolution of the primordial gas density with 
rcdshift assuming z tran = 17 with Az tran = 1.25, and Ztran = 22 with Az t ran = 2.5. The lines show the evolution density 
at constant halo mass. Prior to ztran these follow the expectation from the Press-Schechter mass function. 
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4.5. The Evolution Of The Metal- Enriched Gas Content Of Galaxies 

We compute the density of enriched gas inside halos of mass Mi. In regions of the IGM that have not yet been enriched, 
newly accreted gas is primordial. Conversely in enriched regions of the IGM, newly accreted gas is metal enriched. Prior 
to the reionization of a region, gas can cool inside a halo whose virial temperature is larger than T m i n ~ 2 x 10 4 K. However, 
following reionization of a region, infall is suppressed in halos whose virial temperature is smaller than T ro i on ~ 2.5 x 10 5 K. 
In neutral regions the enriched gas inside halos below M m j n is not photo-evaporated so that, at least initially, some halos 
with masses below M mm will contain enriched gas that has not undergone SF. This gas will move into larger halos during 
subsequent mergers. Using a method similar to that of the preceding section, we can calculate the density of enriched 
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Fig. 2. — Upper panels: The evolution of the density of gas that has not been involved in a starburst, within dark matter halos as a 
function of redshift. The redshift where the IGM became metal enriched was taken to be 2tran = 17, with a variance of Azt ran = 1.25. This 
epoch is marked by the vertical dotted line. The thin curves show the density of gas corresponding to the maximum possible density, i.e. 

" E "M"" h = fbM where /b = O m /Slb- The lines are labeled by the logarithms (base 10) of the corresponding masses. The grey curves 
show the corresponding maximum gas density for the Mj, M m ; n and M re ; on . The thick lines show the evolution of the density of gas that 

has not yet undergone a star-burst. The upper left and upper right panels show the evolution of primordial gas d " p J^° rd , and enriched gas 
jl '™ ch respectively. Lower panels: Collapsed fractions and SF. The lower left panel shows the collapsed fraction of primordial gas (solid line) 
and enriched gas that has not yet undergone a starburst (dashed lines) . For comparison we show the collapsed fractions of dark-matter halos 
above Mj, M m ; n and M re ; on (left to right). The lower- middle and lower-right panels show the SF rates in neutral regions for Pop-Ill stars 
(solid curves) and for Pop-II stars (dashed curves), as well as for Pop-II stars in rcionizcd regions (dashed grey curves). The lower-middle 
and right-hand panels show SF in cases- A and B respectively (see §5). 



gas that remains in these low mass halos following reionization. Note here that we assume enriched gas can accrete into 
halos of all masses in excess of Mj. Conversely, in § 3.1 we argued that the IGM will remain at a temperature of ~ 500K 
following enrichment by a blastwave, implying that accretion of enriched gas is suppressed for mini-halos with masses 
smaller than ~ 5x 10 6 M Q . However this inconsistency will not affect our estimate of Pop-II SF since ~ 5x 10 6 M Q < M m i n . 

We can therefore write down the evolution of the density of enriched gas which we denote dp ^ M i ^ il {note that this density 
refers to a density of gas that has not been in a galaxy at a time when it underwent a burst of SF) . Photo-ionization leads 
to dp , c "', ich = in ionized regions of the IGM. In neutral regions we have 
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otherwise, (13) 

dM \dz 

where & is the Heaviside step function. The first two terms (first line) only contribute in enriched regions of the IGM. 
The first term corresponds to a source of collapsing halos in which metal enriched gas accretes from the IGM. The second 
term corresponds to a source term of gas accumulation due to progenitor halos which are below the Jeans mass. These 
progenitor halos result in accretion of enriched gas from the IGM. The third term describes the movement of enriched gas 
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Fig. 3. — As per Figure 2, but with ^tran = 22 and Aztr. 



2.5. 



from small halos to large halos during mergers. The fourth term describes the loss of gas density in halos of mass M that 
results from the merger of those halos to form larger systems. The enriched gas density in halos smaller than the Jeans 
mass or larger than M cr ;t is zero. The upper right panels of figures 2 and 3 show the evolution of the enriched gas density 
with redshift assuming z tran = 17 with Az tran = 1.25, and z tran = 22 with Az t ran = 2.5. 

4.6. The Collapsed Fraction Of Gas Inside Galaxies 
The fraction of gas in the universe that has collapsed inside galaxies at redshift z but not undergone a starburst is 

F primord (z) = r dM 1 dPp ;;7 d (z) (u) 

Jo /bSWc <*M 



for primordial gas, and 



F..,, i . ] A-.)= i dM 1 dPe ::; h (z). (15) 







for enriched gas. These are plotted in the lower left panels of Figure 2 and 3 assuming Ztran = 17 with Azt ran = 1.25, 
and z tran = 22 with Az tran = 2.5 respectively. For comparison we also show the collapsed fractions (F co \) above Mj, 
M m [ n and M re i on (left to right respectively). We see that prior to the transition redshift, nearly all gas inside collapsed 
systems is primordial. After the transition redshift this collapsed primordial gas is replaced with enriched gas. However, 
the collapsed fraction of primordial gas does not peak until z ~ 14 and z ~ 16 assuming z tra n = 17 and zt ran = 
22 respectively. The enriched gas comes to dominate the collapsed gas fraction at around the same time. Note that 

-^col(-^j) — -fprimord ~t~ -^enrich ~t~ -^col (-^min) ■ 

5. THE STARFORMATION HISTORY 

We model SF to occur following the collapse of baryons into a dark matter halo larger than M m i n , or following a major 
merger with progenitors larger than M m j n . Halos that merge to form a system greater than M m ; n for the first time are 
able to produce Pop-II or Pop-Ill stars. The evolution of gas densities derived in the previous section may be used to 
compute the relative quantities of each. Barkana & Loeb (2000) suggested a model for the SFR based on the merger 
history of halos. In their model, gas that has not previously cooled and undergone an episode of SF in a galaxy forms 
stars with efficiency r\ as it cools following a merger that forms a galaxy larger than M cr j t . In addition, gas in a galaxy of 
mass M > Af m j n undergoes an additional starburst if the galaxy merges with a second galaxy whose mass is larger than 
M/2. Barkana & Loeb (2000) assume that the cold gas in a galaxy has undergone one previous starburst, so that the 
remaining gas mass is reduced by a factor (1 — ry). 

Following Wyithe & Loeb (2003a) we consider two prescriptions for SF which bracket the expected scenario of supernova 
feedback in low mass galaxies. The scenario described above is referred to as Case-A. We also consider a second evolution 
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for the stellar ionizing field, denoted hereafter as case B. Analysis of a large sample of local galaxies shows that the ratio 
e = M*/Mhaio (where M* and Mh a i are the total stellar and dark matter halo masses respectively) scales as e cx M^ a f o 
for < 3 x 10 10 M Q , but is constant for larger stellar masses (Kauffmann et al. 2003). Since SF is thought to be 
regulated by supernova feedback (Dekel & Silk 1986), the important quantity is the depth of the galactic potential well, 
or equivalently the halo circular velocity. Using the stellar mass Tully-Fisher relation of Bell & De Jong (2001), we find 
the threshold circular velocity v± = 176 km s -1 that at z = corresponds to a stellar mass of 3 x 10 10 M Q . In this case 
we define n as the SF efficiency in galaxies with circular velocities larger than v+ . The SF efficiency in smaller galaxies is 
multiplied by a factor e, where e = 1 for M > M^ alo and e = (M/M^ lo ) 2/3 for M < M^ alo . The expressions for the SF 
rate in the following sub-sections have SF efficiencies written as the product en. Note that while in Case-B e shows the 
mass dependence described above, e = 1 for all masses in Case A. 



5.1. The Pop-Ill SF Rate 

We assume that there is a burst of Pop-Ill SF of mass enf^M (where en is the SF efficiency and we assume /b = £!h/f2 TO ) 
whenever there is a merger of a halo with mass M < Mj to form a new halo with mass M\ > Mj. These starbursts 
comprise all of the Pop-Ill SF, since once gas has been involved in a starburst the gas becomes enriched so that a 
subsequent shock due to a major merger results in Pop-II SF. Following reionization of a region of IGM, all gas in halos 
below M m ; n is photo-evaporated over a timescale ~ 0.1-ff -1 where it becomes enriched, and hence subsequent SF there is 
assumed to be Pop-II. The Pop-Ill SF rate density in neutral regions is therefore given by 

f.min(M,M min ) 



- p " m ° rd (z) = / dM dAM^ 

dt JM m i„/2 imax(0,M mb -M) dAMdt 
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ps 
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6(M n 
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dAM 

dn„ B 



(16) 



dAM J / 

The Pop-Ill SF rate in a neutral IGM is shown in the lower middle (Case-A) and lower right (Case-B) panels of figures 2 and 
3 (solid lines) assuming ztran = 17 with Az tran = 1.25, and z tran = 22 with Az tm n = 2.5 respectively. The corresponding 
integrated Pop-Ill SF densities follow from 



Fprimord(^ > ^) 



dt cfTprimord 



(17) 



dz dt 

As noted in Figure 1, the collapsed mass of primordial gas at z tran far exceeds the mass of gas involved in SF up until 
that point. Moreover, we see that in the absence of reionization, Pop-Ill SF can persist long after the enrichment of the 
IGM. Indeed, Figures 2 and 3 show that if enrichment occurred at z tra n = 17 — 22 then Pop-Ill SF would not peak until 
rcdshift z <~ 12, and would not cease until z <~ 4. This result is quite insensitive to the exact value z tran and to the mode 
of SF. 



5.2. The Pop-II SF Rate 

In the model of SF employed, gas that has not previously cooled and undergone an episode of SF in a galaxy forms 
stars with efficiency en as it cools following a merger to form a galaxy larger than M cr ; t . In addition, gas in a galaxy of 
mass M > M m ; n undergoes an additional starburst if the galaxy merges with a second galaxy whose mass is larger than 
M/2. This second star-burst is assumed to be always Pop-II. 

We may use our calculation of the evolution of dp eni i c h/dM to calculate the Pop-II SF rate in neutral regions. 
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dz dM 



X )M + 8(AM — M mill )AM] (18) 

The first term is for SF where gas crossing M m ; n undergoes a starburst. The second term corresponds to the secondary 
starbursts (in gas that has already undergone a burst of SF) that follow major mergers. Note that while the first term 
does not contribute in regions of the IGM that are not yet enriched, this second term can contribute at all regions of IGM 
that have not been reionized. 

Finally, we calculate the Pop-II SF rate in regions of the IGM that have been reionized. 
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dt dn ps 
M dz dM 



l M /2 dAMdt 

xoj(l - er?) [6(M - M min )/ b M + 6(AM - M min )/ b AM] . (19) 

Here the 1st term is for SF where gas in the halo undergoes a starburst when the halo crosses M re i on . The 2nd term 
corresponds to the additional starbursts (in gas within halos above M m ; n that has already undergone a burst of SF) that 
follow major mergers. 

The Pop-II SF rates in a neutral (dark-dashed lines) and an ionized (grey-dashed lines) IGM are shown in the lower 
right panels of figures 2 and 3 assuming z tran = 17 with Az tran = 1.25, and z tran = 22 with Az tran = 2.5 respectively. As 
noted by Barkana & Loeb (2000), the SF rate in an ionized universe implied by a SF efficiency of r\ = 0.1 agrees with 
observational estimates of ~ 0.1M Q /Mpc 3 /yr at z <~ 5. The integrated Pop-II SF density follows from 

F cnrich (>z) = Jjzff^f^iz). (20) 

6. REIONIZATION OF THE IGM AND THE SF HISTORY 

In the previous section we computed the SF rate as a function of time in neutral and ionized regions of the universe 
for both Pop-II and Pop-Ill stars. However, the average SF history of the universe cannot be computed in isolation 
from its rcionization history, which specifies the fraction of SF that is in neutral and reionized regions as a function of 
redshift. Similarly, the computation of a reionization history requires specification of a SF history. As a result the SF and 
rcionization histories must be computed in a co-dependent way. 

6.1. Modeling the Reionization of Hydrogen 

The simplest estimate of the epoch of reionization is based on the following considerations. Given a co-moving density 
of ionizing photons n 7 in a homogeneous but clumpy medium of comoving density no (where the size of the HII region is 
much larger than the scale length of dumpiness), the evolution of the volume filling factor Q of ionized regions is (Shapiro 
& Giroux 1987; Haiman & Loeb 1997; Madau et al. 1999; Barkana & Loeb 2001) 

dQ 1 dn y C dt 

~T = 3 a B —Qn c — , (21) 

dz no dz a 6 dz 

where «b is the case B recombination coefficient, a = 1/(1 + z) is the scale factor, n e is the comoving electron density, and 
C = (nQ)/(n ) 2 is the clumping factor. This equation describes statistically the transition from a fully neutral universe to 
a fully ionized one, and yields rcionization rcdshifts for hydrogen of around 7-12 for fiducial parameters. However, large 
uncertainties arise in both the source term and in the value of the clumping factor (because more rapid recombinations 
lead to a slower evolution of Q). 

A more realistic description of reionization in a clumpy medium is provided by the model of Miralda-Escude et al. (2000). 
In what follows, we draw primarily from their prescription and refer the reader to the original paper for a detailed 
discussion of its motivations and assumptions. The model assumes that reionization progresses rapidly through islands 
of lower density prior to the overlap of individual cosmological ionized regions. Following overlap, the remaining regions 
of high density are gradually ionized. It is therefore hypothesized that at any time, regions with gas below some critical 
overdensity Aj = Pi/{p) are ionized while regions of higher density arc not. The mass fraction Fm(A;) (or equivalently 
A;) therefore evolves according to the equation 

dF M (Ai) 1 dn 7 R(Ai) dt 

dz n dz a 6 dz 

This equation assumes that all ionizing photons are absorbed shortly after being emitted, so that there is no background 
ionizing field, and no loss of ionizing photons due to redshift. We therefore implicitly assume that the mean free path of 
ionizing photons is much smaller than the Hubble length. This should be valid at redshifts not too much smaller than the 
overlap redshift. 

The integration of equation (22) requires knowledge of Py(A). Miralda-Escude et al. (2000) found that a good fit to 
the volume weighted probability distribution for the density as seen in N-body simulations has the functional form 

(A -2/3_ Co) 2- 



P v (A)dA = Acxp 



A- f3 dA, (23) 



2(2V3) 2 

with S = 7.61/(1 + z) and (i = 2.23, 2.35 and 2.48, and C = 0.558, 0.599 and 0.611 at z = 2, 3 and 4. At z = 6 they 
assume [3 — 2.5, which corresponds to the distribution of densities of an isothermal sphere, and solve for A and C by 
requiring the mass and volume to be normalized to unity. We repeat this procedure to find Py(A) at higher rcdshifts. 
The proportionality of So to the scale factor is expected for the growth of structure in an fi m = 1 universe or at high 
redshift otherwise, and its amplitude should depend on the amplitude of the power-spectrum. The simulations on which 
the distribution in Miralda-Escude et al. (2000) was based assumed Q m = 0.4 in matter, — 0.6 in a cosmological 
constant and a$ — 0.79, close to the values used in this paper. 

Equation (22) provides a good description of the evolution of the ionization fraction following the overlap of individual 
ionized bubbles, because the ionization fronts are exposed to the mean ionizing radiation field. However, prior to overlap, 
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Table 1 

Parameters for models presented in Figures 4-9. 
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0.007,0.01 
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0.008,0.30 
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Figure 5 


Case-A 
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20 


100% 


0.009,0.01 


0.009,0.03 


0.010,0.10 


0.009,0.30 


0.20 0.18 0.14 0.11 


Figure 6 


Case-B 
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1.25 


20 


100% 


0.122,0.03 


0.093,0.10 


0.058,0.30 


0.037,1.0 


0.13 0.09 0.07 0.06 


Figure 7 


Case-B 


22 


2.5 


20 


100% 


0.130,0.03 


0.115,0.10 


0.090,0.30 


0.067,1.0 


0.12 0.08 0.06 0.05 


Figure 8 


Case-A 
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1.25 


10 


10% 


0.015,0.01 


0.015,0.03 


0.015,0.10 


0.012,0.30 


0.21 0.18 0.15 0.12 


Figure 9 


Case-A 


22 


2.5 


10 


100% 


0.018,0.01 


0.018,0.03 


0.021,0.10 


0.022,0.30 


0.21 0.18 0.14 0.11 



the prescription is inadequate, due to the large fluctuations in the intensity of the ionizing radiation. A more accurate 
model to describe the evolution of the ionized volume prior to overlap was suggested by Miralda-Escude et al. (2000). In 
our notation the appropriate equation is 

= \*£ - ob(1 + zfR(A c )n e Q^. (24) 
dz n u dz dz 

or 

dQ 1 dn-y 



dz n°F M (A c ) dz 



a B (l + zYR(A c )n 



3 D / a ^ dt , dF M (A c ) 



dz dz 



F M (A C ) 



(25) 



In this expression, Q is redefined to be the volume filling factor within which all matter at densities below A c has been 
ionized. Within this formalism, the epoch of overlap is precisely defined as the time when Q reaches unity. However, we 
have only a single equation to describe the evolution of two independent quantities Q and Fm ■ The relative growth of 
these depends on the luminosity function and spatial distribution of the sources. The appropriate value of A c is set by 
the mean free path of the ionizing photons. More numerous sources can attain overlap for smaller values of A c . Assuming 
A c to be constant with rcdshift, we find that results do not vary much (less than 10% in the optical depth to electron 
scattering) for values of A c ranging from a few to a few tens. At high redshift, these A c correspond to mean free paths 
comparable to the typical separations between galaxies or quasars. Unless otherwise specified we assume A cr ;t = 10 (which 
lies between the values for galaxies and quasars) throughout the remainder of this paper. 

6.2. Effects of Minihalos On Reionization 

The presence of minihalos effects the progression of reionization, by delaying the buildup of the ionizing radiation 
field within ionized regions. Many of the ionizing photons will be used to photo-evaporate the minihalos rather than 
to progress reionization within IGM. Moreover, the finite time over which this evaporation occurs allows SF to proceed 
in halos crossing M m ; n for a short time following the reionization of a region of IGM. In the presence of an intense 
ionizing background radiation field, a large fraction of the gas in halos with virial temperatures below T v j r ~ 10 4 K is 
photo-evaporated (Barkana & Loeb 2002; Iliev, Shapiro, & Raga, 2005; Ciardi et al., 2005) on the halo dynamical time 
(<~ O.IH^ 1 ). Two additions should therefore be made to the above description arising from the effect of minihalos. First, 
we need to take into account Pop-Ill SF within ionized regions due to halos crossing the H-cooling threshold which 
have not yet been photo-evaporated. A simple way to approximate the effect of this additional SF is to reduce the cold 
gas in minihalos, and therefore the resulting SF rate corresponding to crossing the H-cooling threshold by a factor of 
0.5(Q(t) — Q(t — dt))/Q{t) at time t. Here dt ~ O.lif -1 and the pre-factor of 0.5 is approximately the average remaining 
neutral gas fraction in all minihalos which become exposed to ionizing photons between dt and zero years prior to a 
rcdshift z, while the ratio term is just the fraction of minihalos that become exposed to ionizing radiation within the 
preceding dt years. 

Second, while the above consideration enhances the SF in ionized IGM, the minihalos will delay further reionization for 
a period lasting for the photo-evaporation timescale through screening of the ionizing sources by minihalos (Barkana & 
Loeb 2002). We therefore need to estimate the fraction of photons that are absorbed by minihalos rather than contributing 
to the reionization of the IGM. This estimate may be made as follows. An ionizing photon will travel a distance A gas 
through ionized regions of IGM before encountering neutral IGM. On the other hand, an ionizing photon will travel a 
different distance Ahaio before encountering a minihalo that contains cold gas. The fraction of ionizing photons produced 
in ionized regions that reionize neutral IGM rather than intercepting minihalos prior to their evaporation is therefore 
approximated by 

/mini = cxp{-A gas /A nalo }. (26) 

The length A gas may be easily computed within the formalism described. An ionizing photon will typically travel a 
distance A gas = 60kms^ 1 i? _1 [l — Fy]~ 2 / 3 , where Fy is the fraction of the IGM at overdensities below A c (Miralda-Escude 
et al. 2000). On the other hand, to compute the cross-section for interception of a minihalo that contains cold gas we need 
to find the fraction of gas in minihalos at a redshift z that was accreted to within progenitor minihalos within the previous 
interval O.lif _1 . This quantity may be computed from the extended Press-Schechter formalism (Lacey & Cole 1993). For 
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Fig. 4. — Reionization and SF histories for the case where the redshift at which the IGM became metal enriched is taken to be Ztran = 17 
with a width of Azt ran = 1-25 and the metal mixing efficiency is / m i x = 10%. The epoch ztran is marked by the vertical dotted line. The 
top-left panel shows the differential SF rates for Pop-Ill (solid lines) and Pop-II stars (dashed curves). In the top-right panel, we plot the 
cumulative number of ionizing photons per baryon that have escaped into the IGM by redshift z for Pop-Ill (solid curves) and Pop-II (dashed 
curves) stars. The dotted curve marks number of ionizing photons produced by Pop-Ill stars at ztran- The middle left panel shows the 
contributions to the mean metallicity of the IGM from Pop-Ill (solid curves) and Pop-II (dashed curves) SF. The right hand panel shows 
the reionization history as the evolution of the volume averaged (solid line) and mass averaged ionization fractions (dotted line). The lines 
of different thickness show histories with different combinations of (/enrich > /primord) respectively. In each panel the thicker lines correspond 
to larger values of /primord- The lower three panels represent results for series of reionization histories. The value of /enrich that results in 
a final overlap at z = 6.5 (or as near as possible at z > 6.5) is plotted as a function of / pr i mor d in the lower middle panel. The lower left 
panel presents grey-scale and contours showing the level of overlap as a function of redshift and /primord- The grey-scale shows the level of 
ionization (with black representing high ionization). The dotted, dashed and solid lines correspond to Q = 0.2, 0.5 and .999 respectively. In 
the lower right panel we plot the optical depth t cs as a function of / pr imord (g re y line). The dark line in this panel shows the ratio of the 
amount of Pop-Ill SF already completed at ztran to the total Pop-Ill SF. 



a halo of mass M at redshift z, the fraction of the halo mass which by some higher redshift z-i had already accumulated 
in halos is 

r im/D(z 2 ) - i.m/D{z) 



F g (z, Z2) = erfc 



y/2[S(M } ) - S(M)} 



(27) 



where D(z) is the linear growth factor at redshift z, S(M) is the variance on mass scale M (defined using the linearly 
extrapolated power spectrum at z — 0), and Mj{z%) is the Jeans mass at Zi. Given the virial radius for a halo R v i r (M) of 
mass M (Barkana & Loeb 2001) and the Press-Schechter (1974) mass function (dn/dM) for the number density of halos 
we may then find the length Ahaio 



Ahalo — 



M., 



F g (z,z + 5z)} 



(28) 



where Sz = OAH-^dt/dz)" 1 . 
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6.3. Co-Evolution of Reionization, Star Formation and Metal Enrichment 

We next describe the co-evolution of the reionization and SF histories. Having specified a transition redshift, we 

and Pop-II 

dSE.^ TV,, 

dt 

dn y _ dn 7iquasars dt M© 



computed in the previous section the Pop-Ill ( drp ™° rd ) and Pop-II ( <ff '^ teh ) SF rates in neutral regions, as well as the 
Pop-II SF rates in ionized regions ( ' jrcn "° h ''° n ). The ionizing photon production rate then follows from 



dz dz dz ra p 



(1 — Qhoat) ( ^enrich /enrich 3jT^ ( Z ) + -^primord/esc P ^™°' d (z) ) (20) 



r ( f at ^Fenrxchjion / \ . 1 Qhcat Qhcat(^ "F (5z) ^Fp r i moro ! . . 

"T ^Thcat/mini I /enrich^ *cnrich T7 v^/ ~r 77 7^ 'primord/esc T7 V^J 

\ «E ^ tt/hcat «t 

where m p is the proton mass, -/Vp r i morc i and -/Ve nr i c h are the ionizations per baryon, and / pr imord and /enrich are the escape 
fractions for Pop-Ill and Pop-II stars respectively. The quantity Qheat is the maximal value that Q has attained at any 
time in the reionization history. A region that has been reionized may recombine in the absence of an ionizing radiation 
field, but the timescale for cooling is longer than the recombination timescale for > 10. The inclusion of Qhcat ensures 
that the SF rate in regions that have recombined remains the same as for an ionized region should that region return to 
a neutral state. The term d "" ,, ^ ae """ 3 accounts for the contribution of quasars to the ionizing flux. Quasars are included 
in the model as described in Wyithe & Loeb (2003a), with an updated quasar luminosity function model (Wyithe & 
Loeb 2003b). If treated as the sole ionizing sources, we find that quasars reionize cosmic hydrogen at around z ~ 5. 
Wyithe & Loeb (2003c) found that helium became reionized by quasars at z ~ 3.5, in agreement with observation (e.g. 
Theuns et al. 2002). The coupled evolution of equations (22), (25) and (29) yield the SF and reionization histories. 

We may also estimate the mean level of metal enrichment in the IGM. This quantity is important because massive 
Pop-Ill stars are expected to form only out of gas that is not enriched to a metallicity above Z CI - lt ~ 10 Z Q (Bromm et 



14 




al. 2001a). We estimate the mean metallicity at a redshift z as 

Y f°° dF 

ZlGM(z) = /mix / rfz'-TT' ( 30 ) 

w SN J z dz' 

where wsn is the stellar mass that must be produced per supernova, Y is the mass of metals produced per supernova, and 
/ m i x is the fraction of metals that escape from galaxies into the IGM in supernova-driven outflows (e.g. Madau, Ferrara 
& Rees 2001; Thacker, Scannapieco & Davis 2002). The value of / m ix may be large because the early dwarf galaxies had 
shallow potential wells. For a massive Pop-Ill IMF, Furlanetto & Loeb (2003) estimate wsn ~ 462M Q and Y ~ 70.5M Q . 
We note that this number is very sensitive to the stellar mass function, because stars with masses 140 < M < 260M Q are 
expected to eject their metals while more massive stars collapse to form black holes (Heger & Woosley 2002). 

Different reionization histories lead to different column depths of free electrons along a line-of-sight to recombination. 
An important probe of the reionization history is therefore provided by the observation by the WMAP satellite of a large 
optical depth to electron scattering of CMB photons. The optical depth to electron scattering r es , depends on the mass 
filling factor (QFm) and is evaluated according to 

Z" 1000 cdt 

r cs - / dz— a T QF M n° B (l + z) 3 , (31) 
Jo dz 

where ctt = 6.652 x 10~ 25 cm 2 is the Thomson cross-section. 

6.4. Results: The Reionization History of Hydrogen 

Figures 4-9 show several examples of the reionization history, which explore a wide range of the various inputs to the 
model. These inputs are summarized in Table 1. In each case we present 4 histories in the upper set of panels. The input 
values of /primord, together with the value of /enrich that leads to a final reionization at z = 6.5 are also listed in table 1. 
The thicker lines correspond to histories with larger values of / pr imord- Note here that the escape fraction is degenerate 
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with the SF efficiency (for 77 1) . We have chosen r\ = 10% in presenting our SF histories, however the values of escape 
fraction can be adjusted in proportion for different choices of rj. 

The total SF rates in Pop-II (solid lines) and Pop-Ill (dashed lines) stars are plotted in the upper left panels of 
Figures 4-9. Feedback from reionization quenches SF. However, our histories show that the additional Pop-Ill SF remains 
significant at z < z tran in all cases, and peaks at z < z tran in many cases. Pop-II stars begin to dominate the global SF 
rate only around the time of the peak in the Pop-Ill SF rate. We also find that Pop-Ill SF at redshifts below Ztran forms 
the dominant contribution to the total integrated Pop-Ill SF rate in all models. This may be seen through plots of the 
integrated ionizing photon production as a function of redshift (upper right panels of Figures 4-9). The horizontal dotted 
lines show the level of ionizing photon production at z = 2tran- The fraction of ionizing photon production by Pop-Ill at 
z < ztran is larger than 90% in nearly all cases, and can be as large as 99.9%. Furthermore, Pop-Ill SF can dominate the 
total ionizing photon production at z > 6.5 for models with large r cs . We note that the comparison of SF in an ionized 
and in a neutral IGM (see Figures 2-3) suggests that reionization of the IGM should be accompanied by a dip in the SF 
rate as seen in the calculations of Barkana & Loeb (2000). However, the SF histories that are computed in concert with 
a reionization history do not show this dip. The reason is that an extended reionization history washes out this feature 
on average, though it may still be observed in isolated regions. 

For each of the reionization histories in Figures 4-9 we plot ZiGM(ztran) /Zq (middle left panels). We have chosen / m i x 
(see values in Table 1) so as to obtain Z/Zq ~ 10~ 4 ; close to the threshold deduced by Bromm et al. (2001a) at z tran . Note 
that Z/Zq is degenerate with 77 in addition to / m i x - However, this simple estimate demonstrates that metal enrichment of 
the IGM by Pop-Ill stars is consistent with our assumed z tran (see also Mackey, Bromm & Hernquist 2002). The possible 
exception is for Case-B SF with z tran = 22. In this example the integrated SF rate by Zt ran is not sufficient to enrich the 
IGM. Finally, the parameter A c is unspecified a-priori. We have chosen A c = 10, with the exception of Figures 8 and 
9. Comparison of Figures 8 and 9 with Figures 4 and 5 demonstrate that while the details of the reionization history 
depend on the value assumed for A c , the qualitative nature of the predicted histories and the values of r es are similar. 
We therefore find that our conclusions are not sensitive to the choice of A c . 
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We find that either Pop-Ill stars substantially reionize the universe, or that the fraction of Pop-Ill SF which takes 
place following the transition redshift is large. Either way the conclusion is that Pop-Ill stars played a major role in the 
reionization history of the universe. In particular if the universe did not become substantially reionized prior to Ztran, 
then the large majority of Pop-Ill SF takes place at redshifts below z tran , in contrast to previous modeling. Each of the 
histories shown in Figures 4-9 was chosen so that overlap occurs at z ~ 6.5 as suggested by the spectra of high redshift 
SDSS quasars. However, different reionization histories result in different densities of electrons as a function of redshift. 
The values of r os for the various models are listed in the middle right panels of Figures 4-9 and in Table 1. WMAP 
suggested r os > 0.1. We see that only a small value of ry/primord > 10~ 3 for Case-A, or ry/ pr imord ^ 10~ 2 for Case-B 
is required to obtain this large optical depth, even where the enrichment of the IGM occurred at a redshift as high as 
ztran = 22. Indeed the conclusion that Pop-Ill stars contributed significantly to reionization is quite insensitive to the 
transition redshift for ?7/ pr imord ^ 1CP 3 , so long as Pop-Ill stars had a top-heavy IMF. 

The lower three panels of Figures 4-9 represent results for a series of reionization histories. For each value of the 
primordial escape fraction (/primord), we have found the value of /enrich that results in a final overlap at z = 6.5 (or as 
near as possible at z > 6.5), and which maintains the reionization until the present day. In the lower left panel we present 
grey-scale and contours showing the level of overlap as a function of redshift and / pr imord- The grey-scale shows the level 
of ionization (with black representing high ionization). The dotted, dashed and solid lines correspond to Q — 0.2, 0.5 and 
.999 respectively. The value of /enrich corresponding to these histories is plotted as a function of / pr imord in the middle 
panel. In the right panel we plot the optical depth r es as a function of / pr imord (grey line) . For Case-A, the optical depth 
has a maximum near r cs ~ 0.2, but the WMAP result is consistent with all values of / pr imord- F° r Case-B, the optical 
depth has a maximum near r cs ~ 0.1. Only larger values of / pr imord are consistent with the WMAP result for Case-B. 
The dark line in the lower-right panel is the ratio of the amount of Pop-Ill SF already completed at ztran to the total 
Pop-Ill SF. This ratio is small (~ 0.1) for small values of r os . In Case-A, we find that the ratio can be as small as 0.01 
where #tran = 22 and 0.1 where z tran — 17. The ratio is only of order unity if / pr imord is also of order unity. In Case-B, 
we find that the ratio can be even smaller, with values a factor of 3 lower. The ratio does note become of order unity for 
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any value of / pr imord in Case-B. 

Interestingly in Case-A a range of histories show two epochs of reionization. However, contrary to earlier works, the 
tail end of the first reionization has not been set by an abrupt transition redshift z tran . Rather the end of the first 
reionization can come either before, or after Zt ran depending on the sum of the various astrophysics included in the model. 
A particularly interesting history is seen where Case-A SF is combined with a large escape fraction for Pop-Ill ionizing 
photons and a large transition redshift (see Figure 5). In this example the reionization history shows three peaks of 
reionization. Thus, the claim made earlier by Furlanetto & Loeb (2005) that a temporally smooth metallicity enrichment 
of the IGM gives rise to only monotonic reionization process is not borne out in our more detailed calculations in general. 

Clearly, the information embedded in the reionization history is very rich. It is useful to extract some essential properties 
with regard to the CMB experiments such as WMAP, since the three-year WMAP results are expected to be made public 
shortly. The thin black curves in Figures 6 and 7 show histories with r cs = 0.06 and 0.05 respectively. These histories 
follow from a very small Pop-Ill contribution to reionization, suggesting that without Pop-Ill contribution we expect 
Tes < 0-05 — 0.06. At the other extreme, when we maximize the contribution from Pop-Ill stars to reionization by using 
a very high (perhaps implausible) escape fraction (relative to that of Pop-II stars), a maximum value of r cs = 0.21 is 
possible. This suggests that the optical depth reported by the three-year WMAP results should not exceed this value 
(which is ler above the central value of the first year WMAP result), otherwise our understanding of structure formation at 
high redshift would require a dramatic revision. In a more realistic picture where the product of star formation efficiency 
and escape fraction for Pop-Ill stars is comparable to that for Pop-II stars, we find that the range t os = 0.09 — 0.12 is 
most likely, consistent with our previous calculations (Cen 2003a, b). 

The three different outcomes described above imply qualitatively different histories. In the small r es < 0.05 — 0.06 case 
without significant Pop-Ill contribution, reionization is expected to be rapid with the neutral fraction quickly rising above 
50% at z ~ 8. In the large r cs = 0.012 — 0.21 case with a large Pop-Ill contribution, the reionization process could be 
complex and frequently displays multiple peaks. Under seemingly more reasonable assumptions regarding the efficiency 
of Pop-Ill stars, there appears to be a reionization plateau at z — 7 — 12 where the mean neutral fraction holds in a 
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relatively narrow range of 10% — 30%. Independent constraints of the reionization history are highly desirable. Surveys of 
Lya emitters (e.g., Malhotra & Rhoads 2005) and measurements of QSO Stromgrcn sphere sizes (e.g., Mesinger, Haiman, 
& Cen 2004; Wyithe, Loeb & Carilli 2005; Fan et al 2005) may provide powerful constraints on the evolution of the 
ionized fraction close to the end of the reionization epoch. In addition future CMB experiments such as Planck survey 
might be able to distinguish (e.g., Holder et al. 2003) between the various histories and if so, it should yield some critical 
information on SF processes in the high redshift universe. 

7. CONCLUSION 

We show that the collapsed fraction of primordial gas that is trapped in virializcd minihalos at the redshift (z tran ) when 
the IGM becomes enriched with metals is larger than the collapsed fraction of gas that has already been involved in SF 
by a factor of a few to a few tens. We argue that this virialized gas will not be as easily enriched by super-galactic winds 
as the general IGM, because the strong self-gravity created by dark matter and significant ovcrdensity with respect to the 
background gas help stabilize the collapsed gas against large disruptions by sweeping Shockwaves. Moreover, enriched gas 
is expected to have been pre- heated during the process of enrichment, preventing contamination of virialized primordial 
gas during hierarchical growth of mini-halos. Hence, collapsed gas in virialized minihalos may largely remain metal-free, 
until those halos merge to form larger systems of sufficient mass to initiate SF. 

We have followed the evolution of virialized primordial gas through hierarchical mergers of their host dark matter 
halos. This enables the calculation of the Pop-Ill SF rate arising from the virialized primordial gas when it crosses the 
hydrogen cooling threshold. Our model for the reionization of the IGM computes the co-evolution of the reionization 
and SF histories. In addition to following the gas-content of halos below the minimum mass for SF our model addresses 
several astrophysical phenomena. These include feedback from reionization in the form of suppression of galaxy formation 
in a photo-ionized IGM, as well as the photo-evaporation of minihalos (both the screening of ionizing sources and the 
additional low mass SF in ionized regions prior to photo-evaporation). We also allow the IGM to be gradually enriched 
by metals rather than an abrupt transition. 

We find that the Pop-Ill SF rate peaks substantially later (as low as redshift z ~ 10) in cosmic time than the epoch 
at which the IGM became metal enriched, and may extend to z ~ 6. Moreover, we find that the total integrated Pop-Ill 
SF is larger, by a factor of a few to a few tens, than the integrated Pop-Ill SF at the average redshift of IGM metal 
enrichment. This prolonged epoch of Pop-Ill SF has a pronounced effect on the reionization history of cosmic hydrogen. 
Despite the varied feedback mechanisms included in our model we find that the prolonged epoch of Pop-Ill SF leads to 
complex reionization histories which can be multi-peaked under a range of plausible scenarios. 

Observational measurements of the detailed reionization history would provide critical information on the formation 
process of first stars. Specifically, assuming the final overlap redshift to be fixed at z = 6.5 (as implied by observations of 
high redshift quasars), the contribution of Pop-Ill stars to reionization can be quantified and will be tested by three-year 
WMAP results: (1) if Pop-Ill stars do not contribute to reionization, t os < 0.05 — 0.06 and a rapid reionization at z <~ 6 
is expected, with the mean neutral fraction quickly exceeding 50% at z ~ 8; (2) if the product of star formation efficiency 
and escape fraction for Pop-Ill stars is significantly larger than for Pop-II stars, then a maximum r cs = 0.21 is achievable; 
(3) in a (perhaps) more plausible scenario where the product of star formation efficiency and escape fraction for Pop-Ill 
stars is comparable to that for Pop-II stars, t cs = 0.09 — 0.12 would be observed [consistent with our previous calculations 
(Cen 2003a, b)]. In this case the histories would be characterized by an extended ionization plateau from z = 7— 12 where 
the mean neutral fraction stays in a narrow range of 0.1 — 0.3. These results are rather insensitive to the redshift where 
the IGM becomes enriched with metals. 

We thank Jerry Ostrikcr for helpful discussions. This research is supported in part by the Australian Research Council, 
and by grants AST-0407176 and NAG5-13381. 
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